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1 Introduction 

During the last year, the LHC has made major progress in measuring processes with top 
quarks. As a result, the total production cross-section is now known with few percent accuracy 
at both 7 and 8 TeV [1-3] . The precise measurement of the top pair production cross-section 
has allowed high-precision extraction of the strong coupling constant [4]. At the same time, 
the mass of the top quark has been measured [5, 6] with precision matching the one from the 
Tevatron [7]. 

Such a level of agreement between measurements spanning different colliders, collider 
energies and final states unambiguously signifies the commencement of the high-precision 
top quark measurements phase. Given the close relationship between Higgs and top physics 
[8, 9], entering this high-precision phase is particularly significant also in the context of the 
discovery of a Higgs-like particle [10, 11] at the LHC. 

Equally impressive are the theoretical top physics developments of the recent past. During 
the last couple of years a number of calculations with NLO accuracy were performed, that 
accounted for the decay of the top quarks and even off-shell effects [12-20]. Predictions for 
the total inclusive cross-section [21-29] beyond NLO [30-32] were, until recently, exclusively 
based on using NNLL soft gluon resummation [22, 33, 34] and the threshold approximation 
[35] of the partonic cross-section. Very recently also the high-energy limit of the cross-section 
was incorporated in Ref. [36]. 
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As was demonstrated in Ref. [29], see also the discussion in [27, 37], predictions based 
on soft-gluon resummation alone show only modest improvement over the NLO result, since 
effects that are subleading in the soft limit can be numerically as significant. 

The first step towards top pair production in NNLO QCD was undertaken in Ref. [38], 
where the dominant correction from the qq partonic reaction was computed. The remaining, 
numerically subdominant contribution from this partonic reaction was presented in Ref. [39], 
together with the NNLO corrections from the qq, qq' and qq' initiated reactions. 

The results of [38] demonstrate the importance of a complete NNLO calculation and the 
role it plays in reducing the theoretical uncertainty. Motivated by this observation, in this 
work we compute the NNLO correction to the reaction qg — t- tt + X. Our aim is to verify 
the effect of this nominally subdominant reaction and check the quality of the approximation 
[36] derived from the high-energy limit of the qg partonic cross-section. 

The paper is organized as follows: in section 2 we introduce our notation. In section 3 we 
work out the subtraction of the initial state collinear singularities and the evaluation of scale 
dependent terms. The NNLO parton level result for the reaction qg ^ tt + X is presented 
in section 4. In section 5 we discuss the properies of the new NNLO result at parton and 
hadron levels and compare with existing approximations in the literature. In section 6 we 
update our "best" [29] LHC prediction. 



2 Notation 

We follow the notation established in Refs. [38, 39]. At leading power, the total inclusive top 
pair production cross-section factorizes 



C^tot 



d/3 i/3,fi') aij i/3,m',fi') + O(Aqcd) 



(2.1) 



The indices i,j run over all possible initial state partons; /3max = a/1 — ^ni^ / S; \fS is the 
cm. energy of the hadron collider and j3 = y/1 — p, with p = Am? /s, is the relative velocity 
of the final state top quarks with pole mass m and partonic cm. energy ^/s. 
The function $ in Eq. (2.1) is the partonic flux 



1 - /32 



1 /^max ,.2 



1-^2 



expressed through the usual partonic luminosity 



£ij{x,p ) = x{fi(g) fj){x,p ) = X dy dz6{x - yz)fi{y)fj{z) 



(2.2) 



(2.3) 



As usual, are the renormalization and factorization scales. Setting pp = pfi = p, 

the NNLO partonic cross-section can be expanded through NNLO as 



{I3,m^p^) 



dlf + as 



+ 4 



(2) , r ^(2,1) , 7-2 (2,2) 



■ (2.4) 
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In the above equation L = In (/x^/m^), as is the MS couphng renormahzed with Nl = 5 
active flavors at scale fi^ and a^J^'"^^^ are functions only of /3. 

All partonic cross-sections are known exactly through NLO [30-32]. The scaling functions 
0"!^'^^ and (^If^^ can be computed from cr^j\ see section 3. The dependence on fiR / /^f can 
be trivially restored in Eq. (2.4) by re-expressing asifJ-p) in powers of asip-R)', see for example 
Ref. [21]. 

The reactions ij tt + X were computed for i,j = {qq, qq, qq' , qq') through NNLO in 
Refs. [38, 39]. In this paper we compute the NNLO correction to the reaction qg — >• tt+X. The 
only currently unknown contribution to tt production at NNLO is the gg initiated reaction, 
which will be the subject of a future publication. 



3 Collinear factorization and scale dependence 



We follow the setup and notation described in Ref. [39] and denote the collinearly unrenor- 
malized partonic cross-sections as a^\e,p). Then, introducing the functions s-"^ and s-"^ 



defined as s^[^\e,p) = a^^' {e, p)/p and s'f!:\p) = af-\p)/p, the MS-subtracted g'^-initiated 



(n) 

cross-section Sqg reads through NNLO: 



(n) 
3 • • 



5(1) 
9ff 



5(2) 
99 



i'qg ^ 




+ 



2e2 



2e2 1^39 



4°) ^ Pi? + 4°^ ^ P^'J 



19 



1 

Ye 



5(0) 
■^99 



p(i) + pW 

^gq ^ ^qq ^ qg 



(3.1) 



C3 Pi? Pi°) + ^] ® Pi°) ® Pi°) + ® Pi°) ^3 Pi°) + 4°^ ® Pi°) ^3 Pi°) 



99 



99 



99 



99 



99 



99 



99 



^99 



99 



99 



+ 



1 / 1 




Pq^^ + S^qg ® Pgg' + Sgg ^ Pqq' + S^gg ® Pgq 



m j_ 5(1) « p{o) , 5(1) 



(0) 



e V27r 



(3.2) 



with /3o = UCa/G - Nl/3. 

The integral convolutions in Eq. (3.2) are performed numerically, over a set of 80 points 
in the interval /? G (0, 1). The only non-trivial step in this evaluation is the derivation of the 
partonic cross-section s^J through order 0{e). To derive it, we follow the approach of Ref. [32] 
which allows one to derive analytical results for the required partonic cross-sections. The order 



0{e) terms of s^^^ and Sq^g can be easily computed this way and expressed in terms of standard 
harmonic polylogarithms (HPL) [40]. As can be anticipated from the findings of Ref. [32], 
however, the calculation of s^g through order 0{e) introduces a number of new functions that 
go beyond the class of HPL's. In particular, some functions are represented as two dimensional 
integrals. From a numerical point of view, this is problematic since it significantly reduces 
the speed of the numerical integrations in Eq. (3.2). To deal with the loss of speed, we 
have resorted to interpolation techniques, which limits the appeal (and usefulness) of an 



s(i) 
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intrinsically analytic approach. Based on our experience, we conclude that such an approach 
for the computation of the collinear factorization contributions is suboptimal. 

(2 1) (2 2) 

The evaluation of the scale dependent functions Cqg and (iqg is rather straightforward, 



see [39] for details. In terms of the functions s, 
1 



uiuiio uqg aiiu. uqn la lauiic 

(n(,m))/ N _ (n(,rn))f , 



,(2,2) 
19 



+ 



2(27r)2 



get: 



+ 



p(0) 
91 



p(0) 

99 



(0) 



.(2,1) 
99 



1 



(2^)2 

1 

2^ L 



.(0) 
99 



99 ^ ^99 



(0) 



5(1) 
99 



p(0) 

99 



p(0), 
99 

^ 99 



5(1) 
99 



99 



p(0) 

99 



99 



99 



99 



'99 



(3.3) 



Eq. (3.3) agrees with Ref. [21]. The convolutions appearing in Eq. (3.3) are computed 
numerically. We have checked that the fits implemented in the program Hathor [41] agree 

with our own numerical calculation of Eq. (3.3) to a very high precision. Given this level of 

(2 1) 

agreement, instead of producing new fits, we have implemented the analytical fits for aqg 

(2 2) 

and Gqg from Ref. [41] in our program Top++ (ver 1.4) [42]. 



4 Parton level results 

(2) 

For the calculation of the collinearly unrenormalized partonic cross-section dqg we follow 
the approach already used in Refs. [38, 39]. The correction due to double real radiation is 
computed following Refs. [43, 44]. ^ For the real- virtual correction we use the counter-terms 
from Refs. [48-51]. For the evaluation of the required one-loop five-point amplitude we use a 
code from the calculation oi pp — )• tt + jet at NLO [52]. 

As in Refs. [38, 39], Eq. (2.4) is derived in a renormalization scheme where the number 
of active flavors Nf equals the number of light flavors, i.e. Nf = Nl = 5. From a practical 
point of view, the calculation is performed in three steps. In the first step all calculations, 
including UV renormalization, are performed in a standard way by working in conventional 
dimensional regularization (CDR) and considering all fermions as active flavors, i.e. Nf = 
Nl + 1. The renormalization procedure, including the relevant renormalization constants, 
has been described, for example, in [53, 54]. In the second step the heavy flavor is decoupled. 
The decoupling procedure is applied in d = 4 — 2e dimensions, as appropriate, to each of the 
principal contributions to the cross-section: double-real, real- virtual and if present, one- and 
two-loop virtual amplitudes. The decoupling constant can be found, for example, in Ref. [54]. 
The third and flnal step consists of the collinear subtraction described around Eq. (3.3), which 
is performed with all cross-sections (obtained in step two), splitting functions and /3- function 
coefficients evaluated consistently in a scheme with Nf = Nl = 5 active fiavors. 

The result for the NNLO correction to the reaction qg ^ tt + X reads: 

c7g)(/3) = Fo(/3) + Fi(/3)iVi. (4.1) 
^Methods for computing the double real radiation for this process have also been developed in Refs. [45-47]. 
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Figure 1. The functions Fo(/3) and Fi{l3) {the latter multiplied by a factor of 10 for better visibility) 
as defined in Eq. (4.1). Shown are the fits for Fo(/3) (dashed red) (4.3), Fi(/3) (solid blue) (4.2) and 
the discrete computed values, including their numerical errors. 



The full dependence on the number of light flavors Nl in Eq. (4.1) is made explicit. The 
functions Fq^i read: 

Fi = 0.363838/3^ - 1.44391/3^ + 1.1146/3^ - 0.309165/33L/3 + 0.990057/3^L^ + 0.362183^^^ 
+ (0.194867p + 1.57274p2) + 0.0401411pL3 , (4.2) 

Fo = 28.0998/3^ + 24.1753/3^ - 12.3211/3^ - 49.909/3'^ + 11.7853/33l^ + 28.6697/3^L^ 

+ (-1.68957 + 30.6335p2) Lp + (-9.80339/9 - 76.7407/9^) - 3.82993pL^ , (4.3) 

where Lp = ln{p), = ln(/3) and we recall that p = 1 — fi"^. 

The functions Fq^i in Eqs. (4.2,4.3) are fits to the numerically computed partonic cross- 
section. The calculation of the function Fq is done in 80 points in the interval /3 G (0,1). 
The highest computed point is /^so = 0.999. The fit and the computed points, including 
their numerical errors, are shown in fig. 1. Except for the very last point /^SO) the quality 
of the calculation is high, sub-1%. The quality of the fit is also good; it fits the computed 
points within the numerical uncertainties for large and moderate (3. For smaller values of (3 
the quality of the fit is not as high, yet the relative deviation of the fit from the computed 
mid-points is better than 1%. Only in the region of very small /3 the relative distance between 
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the fit and computed central values surpasses 1% but in that range the absolute size of the 
result is negligible which makes such deviation phenomenologically irrelevant. 

The calculation of the function Fi is done in 81 points in the interval /3 € (0, 1). In 
addition to the 80 points used in the computation of Fq, we have added the extra point 
/3 = 0.99999. The fit and all computed points, including their numerical errors, are shown 
in fig. 1. The quality of the fit is high, sub-1%, for /3 > 0.1. In the region of smaller /3 the 
quality of the fit deteriorates, but the absolute difference between the fit and computed points 
is extremely small and also phenomenologically irrelevant. 

(2) 

The most prominent feature of the partonic cross-section aqg is its high-energy behavior 
[30, 55-59]: 

^^'J^u-,x\^^, « Hp) +C0 + 0{p) . (4.4) 
The constant ci has been predicted exactly in Ref. [60], with A^^^-independent numerical value 



ci = -1.689571450230512 . (4.5) 

To improve the high-energy endpoint behavior of the fits (4.2,4.3), we have imposed on 
them the exact ~ ln(p) behavior from (4.4). Then, from the fits (4.2,4.3), we derive an 
estimate of the constant cq appearing in Eq. (4.4). We get the value: 

Co = -9.96 + 0.0345iVi . (4.6) 

Setting Nl = 5 we find that Eq. (4.6) agrees ^ with the numerical estimate of cq derived 
in Ref. [36] with the help of completely independent methods. 

The numerical error on the proportional to Nl term in Eq. (4.6) is likely rather small, 
thanks to our ability to extend the calculation of Fi to (3 as high as (3 = 0.99999 and to the 
fact that the function Fi behaves ~ const at large (5. 

On the other hand, estimating the error on the iV/,-independent part of Eq. (4.6) is 
much harder. The reason for this is that the region below = 0.999 (which is the highest 
computed point for Fq) is still not close enough to the high-energy endpoint to be dominated 
by the high-energy expansion (4.4). Going beyond the highest computed point /Jgo = 0.999 
is currently unfeasible since the computational cost for a single point, located well above the 
point /Sgo, would be comparable to the computational cost for all 80 calculated points. 

Combining the above observations with the fact that the numerical error in this last 
computed point is larger, exceeding 1%, we conclude that the error on the A''£,-independent 
part of Co could possibly be as large as few tens of percent. 



5 Discussion 



5.1 Properties of the parton level result 

The most striking feature of the 0{ag) correction to dqg is the similarity of its shape and size 
to the long known 0{a^g) correction. ^ In fig. 2 we compare the two, including the appropriate 

■^We note that the prediction for the constant Co derived in Ref. [36] contains no explicit Nl dependence. 
^For short, in the rest of this section we refer to these two corrections as NNLO and NLO, respectively. 
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Figure 2. Comparison of tlie NLO and NNLO corrections to tlie partonic cross-section, including the 
relative power of as(rnt) as in Eq. (2.4): ascTq]]' (dashed blue) and a's<Jq'f (solid red). 

relative powers of the strong coupling as{mt) ~ 0.1068, see Eq. (2.4). We observe that the 
main difference between the two curves is in their high-energy behavior, which is more singular 

(2) 

in the case of aqg . The similarity in size and shape between the two consecutive perturbative 
corrections indicates that large perturbative NNLO corrections can be expected. Indeed, if 
it was not for the suppression due to the additional power of as, the NNLO correction could 
have been even more sizable, a feature that might be relevant for the description of lighter 
fermion pair production, like bottom quarks. 

To better assess the phenomenological significance of the similarities and differences be- 
tween the NLO and NNLO corrections, in fig. 3 we plot their product with the partonic 
fluxes for the Tevatron and LHC 8 TeV, see Eq. (2.1) for precise definition. In all cases we 
use MSTW2008nnlo68cl pdf set [61]. The relative powers of as{mt) are also included. 

We observe that the similarity in shape and size between the two corrections is preserved 
at the Tevatron. Therefore, one can anticipate NNLO contribution to the hadron-level cross- 
section (Jtot that is similar in size to the NLO one. On the other hand, at the LHC 8 TeV, we 
observe a dramatic difference between the shapes and sizes of the NLO and NNLO corrections. 
Clearly, at this particular collider energy, one can anticipate very strong cancellation between 
the positive and negative nodes of the NLO result, while the NNLO one stays mostly negative. 
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Figure 3. NLO and NNLO corrections to the partonic cross-section (as in fig. 2) times the flux at 
Tevatron (left) and LHC 8 TcV (right), see Eq. (2.1). 

We also note that the high-energy rise of the NNLO correction is completely screened by the 
flux which vanishes in the same limit. 

To better quantify the differences between the NLO and NNLO corrections, in table 1 
we present their separate contributions to citot- The results on lines 1 and 2 as well as Utot 







Tevatron 


LHC 7 TeV 


LHC 8 TeV 


LHC 14 TeV 


h 


Due to a^qg [pb] 


-0.068 


-0.88 


-0.48 


9.01 


h 


(2) 

Due to Uqg [pb] 


-0.057 


-1.82 


-2.25 


-4.07 


h 


4g^(Hathor;(A + S)/2) [pb] 


0.040 


5.78 


8.11 


27.36 


h 


(/3 - h)/otot [%] 


1.4 


4.9 


4.7 


3.7 



Table 1. Central values for the contributions of a\g and aqj to crtot for the Tevatron and LHC 7,8 
and 14 TeV. Also shown is the corresponding contribution from the program Hathor based on the 
approximation to aqg of Ref. [36]. Line 4 shows the difference between the exact result (this paper) 
and the approximation from Ref. [36], relative to utot- 



on line 4 are computed with version 1.4 of the program Top++ [42] with default settings, 
mt = 173.3 GeV, central scales and MSTW2008nnlo68cl pdf set [61]. The numbers in line 
3 are computed with the program Hathor [41] using the same parameters as above. We 
comment on the results obtained with Hathor in section 5.3. 

Prom the first two rows of table 1 we conclude that the absolute size of the NNLO 
correction can be significantly larger than the NLO one depending on the cm. energy of 
the collider. Comparing to fig. 2, however, we note that the large differences between the 
0(a|) and 0{ag) corrections seen in table 1 do not necessarily indicate a breakdown of the 
perturbative expansion, since they result from large accidental cancellations that strongly 
depend on the collider energy. 
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Next, in fig. 4, we compare the exact result for aqj witli its leading-power high-energy 
approximation (4.4). In complete analogy with the case of fermion pair initiated top pair 
production [39] we observe that the high-energy approximation is justified only very close 
to the high-energy endpoint and is a poor approximation to the exact result outside of this 
narrow range. 



5.2 Properties of the hadron level result 

To better judge the effect from the inclusion of the NNLO correction to the qg reaction, in 
table 2 we give the central values for our best prediction for the Tevatron and LHC 8 TeV in 
the following cases: 

• With, or without, soft gluon resummation: the gg reaction is included in NLO+NNLL 
or in approximate NNLO (defined as in Refs. [29, 38]), while all other reactions are 
included in NNLO+NNLL or NNLO. 

• With, or without, the NNLO corrections to qg ^ ti + X. 
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Besides the central values, in table 2 we also show the average scale uncertainty, defined as 
(scale_|_ + scale_)/2. As elsewhere in this article we use MSTW2008nnlo68cl pdf set [61], 
mt = 173.3 GeV and scale and pdf variations are performed as described in Ref. [29]. 





Resummed; 


Resummed; 


Fixed Order; 


Fixed Order; 




gg-included 


no qg 


gg'-included 


no qg 


Tevatron: central [pb] 


7.010 


7.067 


6.949 


7.006 


Tevatron: aver, scale var. [%] 


±2.6 


±2.7 


±4.4 


±3.7 


LHC 8: central [pb] 


220.4 


222.7 


218.5 


220.8 


LHC 8: aver, scale var. [%] 


±5.3 


±7.3 


±4.5 


±6.5 



Table 2. Central values and average scale variations, the latter defined as (scale_|_ + scalc_)/2, of our 
"best" prediction for the Tevatron and LHC 8 TeV. Numbers are given for the following four cases: 

(2) 

with/ without soft gluon resumniation and with/ without including (Jqg ■ The values of the various 
parameters used in the calculation are specified in the text. 



We observe that the effect on the central value from the inclusion of the NNLO qg 
correction is moderate, and brings down the central value by about 0.8% at the Tevatron and 
by about 1% at the LHC. Such a shift is perfectly consistent with our estimate [38] of the 
theoretical uncertainty at the Tevatron. We also note that, as might be anticipated, the size 
of the shift in the central value is the same independently of the inclusion (or not) of soft 
gluon resummation. 

The effects of the NNLO qg correction on the size of the scale variation is more conse- 
quential. At the Tevatron, the uncertainty in the pure fixed order prediction increases by 
about ±0.7%, while the uncertainty in the soft gluon resummed result is unaffected by the in- 
clusion of the NNLO qg correction. This is consistent with the expectation that the dominant 
source of uncertainty at the Tevatron is already accounted for. The effect on the fixed order 
prediction is also at a level similar to the anticipated [38] NNLO correction in gg ^ tt + X. 

At the LHC, on the other hand, we notice a dramatic ±2% decrease in scale uncertainty 
both with and without including soft gluon resummation. This is a significant improvement 
in the precision of the theoretical prediction at the LHC. Despite this improvement, however, 
it is clear that the unknown genuinely NNLO correction in the gg reaction still dominates the 
uncertainty at the LHC. This is evident, for example, from the fact that the scale variation 
of the resummed result is larger than that of the fixed order result (which we take as a more 
conservative estimate of the theoretical uncertainty [29]). 

Before closing this section we address the question of how the uncertainty in the derived by 
us constant cq (4.6) propagates into phenomenological predictions. As we argued in Ref. [38], 
the most natural way to address this question is to consider the ratio: 



(5.1) 



where: 

^qg{P80)= dP%g{(3)ag\p). (5.2) 

The meaning of the function S^g is as fohows: when its argument is the highest computed 
point Pso (/^so = 0.999 in the case of the function Fq), the function T,qg contains the complete 
contribution to dtot due to the part of the fits (4.2,4.3) that is beyond the highest computed 
point, i.e. from the region where our calculation is not derived but extrapolated. As a 
conservative estimate we take the case of LHC 14 TeV where the partonic fluxes are most 
enhanced in the high-energy region. We find that iJ^g (0.999) ^3x 10"^ which is completely 
negligible. In this sense, the uncertainty on the derived constant cq is of no phenomenological 
significance for top pair productions at the Tevatron and LHC. However, for applications of 
Eqs. (4.2,4.3) to top pair production at future higher energy hadron colliders or for lighter 
quark production, like bottom quarks, it would be advisable to re-assess the smallness of the 
ratio (0.999). 

5.3 Comparison with existing approximations 

In most past studies of top pair production beyond NLO, the qg reaction has received little 
attention, and its NNLO correction has, typically, been neglected. To that end it would be 
interesting to compare the exact result derived in this paper with the only approximation to 

(2) 

(Tqg derived previously [21, 36] and implemented in the program Hathor [41]. 

In fig. 5 we plot the envelope of predictions for Uqg introduced in Ref. [36] (grey band). 
The prediction of that reference is based on matching the threshold term introduced in 
Ref. [21] to the high-energy behavior of the cross-section. The spread of the predictions 
reflects the uncertainty in the prediction of the constant cq as estimated in Ref. [36] . ^ The 
blue curve in fig. 5 denotes the "best" approximation of Ref. [36]. While not explicitly shown 
in fig. 5, the threshold term [21] is essentially identical to the blue line in the region (3 < 0.5. 

It is obvious from fig. 5 that, except in the limit of extremely large /3, the exact result 
for Gqg derived in the present work (red line) has qualitatively different behavior compared 
to the approximation of Ref. [36]. The observed disagreement applies also to the threshold 
term introduced in Ref. [21]. Given the significance of these differences, it is imperative to 
quantify their phenomenological impact. 

(2\ 

In table 1 we present the contribution to cJtot of Ref. [36] 's "best" approximation to Uqg 
(denoted as (A+B)/2 " and corresponding to the blue line in fig. 5). For its computation 
we use the program Hathor [41] with the same parameters as in the rest of this paper. 

We observe that the difference between the exact result and the approximation of Ref. [36] 
is numerically significant at all collider energies. In particular, at the LHC with cm. energies 
of 7 and 8 TeV, the relative difference (with respect to citot) between the approximation of 
Ref. [36] and the exact result can be as large as 5%. Such a shift in the total hadronic 

■'We remind the reader that both our results and the results in Ref. [36] have the same logarithmic behavior 
in the high-energy limit. 
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Figure 5. Comparison of tlie exact partonic cross-section aqg^ (solid red) witli tlie approximation of 
Rcf. [36] (grey band). The central value (dashed blue) is the "best" approximation of Ref. [36]. 



cross-section is very large given that (a) it originates in a subleading channel and (b) it is 
comparable in size to the total theoretical uncertainty at the LHC. We are therefore led to 
the conclusion that such large discrepancy is calling into question the usefulness of the high- 
energy approximation of the heavy flavor production cross-section as a means of describing 
top pair production at hadron colliders. 



6 Phenomenological predictions 

Implementing the 0{ag) correction to aqg (4.1) in version 1.4 of the program Top++ we obtain 
the following "best" predictions for the Tevatron and LHC 8 TeV: 

^NNLO+NNLL(Te,atron) = 7.010 ^_lif,^^Si [scales] t^l^^^J,^], [pdf] , (6.1) 

^(N)NLO.NNLL(LHC3,,^) = 220.4 ^-.^g.^J) [scales] ^^gg [pdf] . (6.2) 

Theoretical prediction for any other LHC cm. energy can be easily obtained with version 1.4 
of the program Top++ by adjusting the collider energy in its default LHC setting. 
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Figure 6. Our "best" prediction for tire LHC as a function of the collider energy. Inner band (light 
green) represents scale uncertainty; total band is the linear sum of scale and pdf uncertainties. Also 
shown are the most precise measurements from CMS and ATLAS [1-3]. 

The numbers above are derived in the following way: the partonic reaction gg ^ ti + X 
is included at NLO+NNLL, as in Ref. [29]. All other partonic channels are now known in full 
NNLO and are therefore included with the exact NNLO results, including NNLL soft gluon 
resummation for the qq ^ tt + X reaction. We use MSTW2008nnlo68cl pdf set [61], and 
scale and pdf variations are performed as described in Ref. [29]. 

We present our predictions for mt = 173.3 GeV. Such value for mt is consistent with the 
current best measurements from the Tevatron [7] (173.18 ± 0.94 GeV), CMS [5] (173.36 ± 
0.38 ±0.91 GeV) and with the ATLAS and CMS top mass combination [6] (173.3 ± 1.4 GeV). 

The measurements [1-3] we compare to, are presented at mt = 172.5 GeV (both for 7 
and 8 TeV). For a consistent comparison, we translate all measurements to mt = 173.3 by 
rescaling them with a common factor of 0.993512 [1]. In principle each measurement should 
be rescaled with its own scaling factor, however such rescaling is available only for Ref. [1]. 
Given the week dependence of the measurements on the value of the top mass, however, any 
inconsistency due to this procedure is at the sub-percent level and is thus inconsequential 
given the size of the experimental and theoretical uncertainties. 

As we discussed in detail in section 5.2, the inclusion of the NNLO correction to the 



qg reaction has notable impact on the scale dependence at the LHC. For this reason, with 
this paper, we update our NLO+NNLL LHC prediction from [29]. In fig. 6 we compare 
our best prediction (6.2) with the most precise measurements from CMS and ATLAS [1-3]. 
We note a very good agreement between theory and data at both 7 and 8 TeV. At 8 TeV 
the total theoretical uncertainty is comparable to the experimental one, while at 7 TeV the 
experimental uncertainty is almost a factor of two smaller than the total theoretical one, 
mostly thanks to significantly reduced systematics. We are hopeful that the inclusion of the 
full NNLO correction in gg ^ tt + X in the near future will further reduce the theoretical 
error. 

We also calculate the ratio of the cross-section evaluated at 8 TeV and 7 TeV. We find 
that the central value of the ratio, and its uncertainty (evaluated as a restricted scale variation 
of the ratio, see [9]) are not significantly different from the numbers reported in Ref. [9]. 

Before concluding this section we note that the large NNLO correction in qg ^ tt + X 
(relative to the NLO correction in the same reaction) could be indicative of this reaction's 
possible relevance to the resolution of the ApB puzzle at the Tevatron [62 , 63] . 

7 Conclusions 

In this paper we calculate the NNLO (i.e. 0{ag)) correction to the total top pair production 
cross-section in the partonic channel qg — )• tt + X. We follow the computational approach 
already used in Refs. [38, 39] and compute the partonic cross-section numerically in 80 points 
on the interval /3 G (0, 1). The numerical precision of the calculation is high, typically below 
1%. For the practical implementation of the result we have derived analytical fits that have 
simple analytical form. 

Our result is consistent with its expected endpoint behavior: it vanishes at threshold 
/3 = and diverges logarithmically in the high-energy limit /3 = 1. By imposing the known 
[60] exact logarithmic behavior in the high-energy limit we extract the constant in the leading 
power term. The value of this constant agrees with a recent prediction in Ref. [36]. While 
the uncertainty spreads in each of the two results are not small, the observed agreement is 
nevertheless an important consistency check on both setups. We have demonstrated that 
the uncertainty on this constant is completely irrelevant phenomenologically for top pair 
production at the Tevatron and LHC. 

The phenomenological impact of the NNLO qg correction is moderate. At the Tevatron 
its only effect is to lower the prediction of Ref. [38] with approximately 0.8% which is well 
within the total theoretical uncertainty. The inclusion of the NNLO qg correction at the LHC 
lowers the cross-section by approximately 1% while at the same time it decreases the scale 
uncertainty by about ±2%. This is a significant improvement in the theoretical prediction, 
which agrees well with the most recent LHC measurements at 7 and 8 TeV. 

At present, the dominant source of theoretical uncertainty at the LHC is the lack of the 
genuinely NNLO correction in the (^f^-initiated reaction. We hope to report results for this 
last missing at NNLO channel in the near future. 
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